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Abstract 

Precision measurements at the Z resonance agree well with the standard model. 
However, there is still a hint of a discrepancy, not so much in Rb by itself (which has 
received a great deal of attention in the past several years) but in the forward-backward 
asymmetry A'p^ together with R},. The two are of course correlated. We explore the 
possibilty that these and other effects are due to the mixing of 6^ and bn with one or 
more heavy quarks. 



Ever since the Z boson was produced as a resonance at the e~e^ colhder LEP at CERN, 
precision measurements of electroweak parameters as well as as of quantum chromodynamics 
(QCD) became available. Certain deviations from the predictions of the standard model have 
been observed in the past, notably the excess in 

R, = ^ (1) 

Y{Z ^ hadrons) 

This had prompted a flood of theoretical speculations regarding the possible existence of new 
physics]^]. At present, however, the experimental data[0| have settled down to a value of i?6 
consistent with an excess of only 1.3cr and it is certainly not an indication of new physics 
by itself. On the other hand, the forward-backward asymmetry J^p^ is now measured to be 
—2.0a" away. [This quantity used to be less accurately measured and it was always less than 
±1.0(7 from the standard- model prediction.] If one takes seriously the two measurements 
together, a possible discrepancy still remains. In this paper we will explore how the mixing 
of 6l and 6/? with one or more heavy quarks would explain the present data. 

In the standard model, using = 175.6 ± 5.5 GeV and assuming itlh = 300 GeV, the 
overall best fit givesU sin^ 9w = 0.23152, for which Rb = 0.21576 and = 0.10308. The 
experimental measurements are0, Q 

Rb = 0.2170 ±0.0009 (±1.38a), (2) 
= 0.0984 ±0.0024 (-1.95(t), (3) 

where the number of a's is the "pull" which is defined as the difference between measurement 
and fit in units of the measurement error. 

Consider the couplings of the b quark to the Z boson: 

gZ^ 

where the subscripts L,R on b refer to the left and right chiral projections (1 =i= 75)72 respec- 
tively. Hence Rb (x g\ + g\, whereas /^^pQ oc {g\ — gj^j/igl + gji). /,From the data, it is clear 
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that a larger gj^ is desirable because that would decrease A^J^^ and increase Rb. Previous 
attempts to increase Rb mostly considered increasing g1. Two specific exceptions]^, 
proposed to increase gj^ and we will discuss them in detail below. 

In the standard model, 

In Fig. 1 we plot Rb versus A'p!'^ as a function of gjiigl) with gligji) fixed at its standard- 
model value. The experimental range is also displayed. For gj^ fixed at {gL)sM, the best fit 
is 

gl = 0.00736 (6) 

for which Rb = 0.2174 and = 0.1015 are obtained. If we let both g1 and gj^ be free 
parameters, then we get the central values of Rb and with 

gl = 0.17586, gl = 0.00994. (7) 

We now discuss how the above two cases, i.e. Eqs. (6) and (7), may be obtained. In 
Ref. [6], a vector doublet of quarks with the conventional charges, i.e. 2/3 and —1/3, is added. 
We call this Model (A) with (Qi, Q2)l,r ~ (3, 2, 1/6) under SU{3)c x SU{2)l x U{1)y- Since 
{Qi,Q2)l transforms in the same way as the known quark doublets, we define it precisely 
as the one that forms an invariant mass with {Qi,Q2)r- Hence the mass matrix linking 
Q2l) with (bji, Q2r) is given by 



Ml 



n 



nib 
rriQb M 



which shows that bfi-Q2R mixing is dominant, and that 6l-Q2L mixing is suppressed by 
mb/mq and is thus negligible ||ll|. We now have 



9l 



1 /II 
sin^ 9w cos^ 6*2 + ( \ — sin^ 9iv 



3 V 2 3 





2 r 


) sin^ 02 





1 1 

- sin^ 9w sin^ 02 

3 2 



n 2 



(9) 



In order to increase g]^ from its standard-model value, it is clear that sin^ 02 must be greater 
than (4/3) sin^ 9w- Hence a rather large mixing with Q2 is required in this model. Numeri- 
cally, to obtain Eq. (6), we need 

sin^ 02 = 0.3260. (10) 

In Ref. [7], a vector doublet of quarks with the unconventional charges —1/3 and —4/3 is 
added. We call this Model (B) with (Qa, Q4)l,r ~ (3, 2, —5/6). The b-Q^ mass matrix is of 
the same form as Eq. (8) because there cannot be a biQaR term for lack of a Higgs triplet. 
In this case. 



9l 



1 





2 r 


) singes 





sin^ Ow + 7: sin^ 6^ 



(11) 



Now we need only a small mixing to obtain Eq. (6), namely 

sin^ 03 = 0.0173. 



(12) 



For comparison against the above two vectorial models, we consider also the addition of 
one mirror family of heavy fermions. The heavy quarks here are right-handed doublets and 
left-handed singlets. We call this Model (C) with (Q5, Qq)r ~ (3, 2, 1/6), Q5L ~ (3, 1, 2/3), 
and Qql ~ (3, 1, —1/3). The b-Qe mass matrix is then 



(13) 



rriQb ruQ 

which allows both bn-Qen and b^-QeL mixings, so that Eq. (7) may be satisfied. However, 
it is a somewhat unnatural solution because mi, and rnq come from the vacuum expectation 
value of the Higgs doublet, whereas m^g and mQb are invariant mass terms. It is thus difficult 
to understand why the latter two masses are not much greater. Using 

2 r 1 1 I n2 



9l 



1 1 

■2 + 3 



^ + ^ sin^ Ow ] cos^ Ogl + ^ sin^ Ow sin^ O^l 



1 1 

■2 + 3 



^ + ^ sin^ 9w + - sin^ O^l 



(14) 



and Eq. (9) with 6*2 replaced by O^r to fit Eq. (7), we find 

sin^ e^L = 0.0189, sin^ Ogr = 0.3537. 



(15) 



In Model (A) and Model (C), large mixing of Br with a heavy quark is required, as shown 
in Eqs. (10) and (15) respectively. This has important implications on the electroweak oblique 
parameters S,T,U or €1,62, e^. In Model (A), assuming that Qi does not mix with t, c, or 
u, we have the following physical doublets: 



Qi\ I Qi 

Q2 / ^ ' y Q2COS92 - b sin 6*2 / ^ 



(16) 



which would contribute to T or ei. In the above, the masses of Qi and Q2 are related by 
mi = 1712 008 02, assuming that << mi, 1712- Let x = sin^ 62, then we find 

167r sm Ow J^w 



where 

ln(l — x) 



F{x) = -2(1 -x)(2 -x) 



1 + 



2x + 3x1 



X 

Note that F(0) = and F(l) = 1 as expected. Also, F{x) > for < x < 1. Taking 
X = 0.3260 as in Eq. (10), we get Fi(x) = 0.141. Let us choose mi = 200 GeV so that the 
decay Qi b + W would not be a significant contribution to the top signal at the Tevatron. 
In that case, 1712 = 244 GeV and Aei = 2.6 x 10^'^ which would take this model far awaylQ 
from the data. Since our purpose is to find out if mixing with heavy quarks would improve 
the overall agreement with data, this numerical result tells us that Model (A) as it stands is 
not the answer. 



In Model (C), the physical doublets are 



Q5 \ / t 

Qa cos 0QR - b sin Oqr J \b cos 6^1 - Qe sin e^L ^ ^ 



(19) 
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but since O^l is small, it can be neglected, and is unrelated to mg. We now find 

3a; 1 



167r sin^ Ow 



w 



■ 2m|m|cos^ mr 
+ ttIq cos uqji 5 5 m — ^ 



2 2 

m5 — rriQ 



(20) 



Hence we can fine-tune m\/m\ to make Aei small. For example, if we let = 200 GeV, 
then the above expression is minimized with rriQ = 273 GeV for which Aei = 0.52 x 10^^. 
This much smaller shift is acceptable. On the other hand, unlike Models (A) and (B) where 
the heavy quarks are doublets in both left and right chiralities. Model (C) has Q^l and Qql 
as singlets, hence the shift in 5* or 63 becomes nonnegligible. Let x = sin^ 9qr and assume 
Mz « m^jUiQ, then we find 



a 



247rsin2^ 



w 



3 - 8x + Sx - In ^ - x{2 - 3a;) In 



mt 



1.8 X 10"^ (21) 



for X = 0.3537, = 200 GeV, and mg = 273 GeV. This shift would already take this 
model far away0 from the data, not to mention that there is also the leptonic contribution 
of 0.44 X 10^^. Hence Model (C) is also not the answer. 

Let us go back to Model (A) and try to reduce Aei of Eq. (17) by allowing Qi to mix 
with t. In that case, we have 

QicosOiR — tsinOiji \ I QicosOil — tsmOn \ I tcos6iL + QisinOiL 



Q2 cos 92R - b sin 



Q2 



where the masses of Qi, Q2, and t are related by 



nil 



= M{cos9il/ cosOifi), m2 = M / cos 92R, mj = M(sin 6'il/ sin 6'i/j) 



(22) 



(23) 



where M is defined as in Eq. (8). After a straightforward calculation, we find 

3a 1 



Aei 



87r sin^ Ow 



+Aim\ In m\ + A2m\ In m\ + Atm\ In ml 



where Cir = cos6'iij. Sir = sin 6*1/2, ete., and 

A _ 2 2 2_i_4_i_('/|2_2 2_2^ 

^1 ~ '^1R^2R ~ ^IL ' ^IR ' y^'^lR '^1R'^2R '^IL, 



2 

mf 



ml 



iR'iR^^2 



ml 



2 ' 
mf 



(24) 



(25) 



A — 2 , / A 2 2,2 2|2\ 
^2 — -52^ -t- I, ^C^l'^2R + ^IR^2R + ^IL) 



A _ i 22 2 /22\ 

^* — Cj^R ^1R^2R ^IL [^IR^IRJ 




2 ' 



(26) 



(27) 



Note that Aiml + A2ml + Atm^ = as expected. We now let mi = 200 GeV and using 
rrit = 175.6 GeV, we have 



For a given value of c^^, we then fix c^^ and hence M (= itliCir/cil) as well as 1712 (= 
M/c2r), assuming of course that = 0.6740 from Eq. (10). We vary c^^ and compute the 
right-hand side of Eq. (24) numerically. We find that it is in fact a monotonically increasing 
function of decreasing c^^. Hence the value obtained earlier for Aei assuming no Qi mixing 
{i.e. c\ii — c\j^ — 1), which was already too far away from the experimental data, cannot be 
reduced, and Model (A) is not saved by additional mixings. 

Now that both Models (A) and (C) are eliminated by the precision electroweak measure- 
ments, we focus our remaining discussion on Model (B). The additional heavy quarks 
and Qi have charges —1/3 and —4/3 respectively. The physical doublets are 



is a singlet. Hence Q,i decays into h + and would have been observed in the top-quark 
search at the Tevatron if its mass is below 200 GeV or so. In fact, the Tevatron top-quark 
events cannot tell h from 6, hence it is even conceivable that was actually discovered 
instead of t. However, Q3 would also have been produced, since ma = m4/ cos ^3 ~ 1.01m4, 
and it would decay into either b + Z or b + H. The nonobservation of the b + Z mode at 
CDF requires the b + H mode to dominate. However, since LEP data already require the 




(28) 




(29) 



whereas 



bn cos 63 + QsR sin 63 



(30) 
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Higgs scalar mass to be greater than about 90 GeV, the b + H mode cannot dominate and 
this exotic possibihty is ruled out. We conclude that Q3 and Q4 are hitherto undiscovered 
and must be heavier than about 200 GeV. Note that this is beyond the reach of LEP for 
producing bQ-^ + bQ^- 

We have so far assumed that Q3 mixes only with b, but of course it could also mix with 
s and d. In that case, the state b in Eqs. (29) and (30) should be considered as a linear 
combination of b, s, and d, but dominated by b. As a result, there would be fiavor-nondiagonal 
couplings of the Z to ds + sd, db+bd, and sb+bs, and in addition, the charged- current mixing 
matrix mediated by W would lose its unitarity. These couplings are presumably very small, 
but they could affect the standard- model phenomenology regarding the — K^, — B^, 
and B^ — B^ systems. Furthermore, the indirect effect of fiavor-nondiagonal neutral currents 
in the right-handed sector could result in the failure of the Standard Model to describe all 
data, especially the precision measurements to be obtained with the upcoming B factories, 
at KEK and at SLAG. For example, the exotic contributions to the decay B X^'y in vector 
quark models [including our Model (B)] have been analyzed recently The dominant extra 
contribution is from the violation of unitarity in the charged- current mixing matrix. It 
gives rise to nontrivial constraints on the off-diagonal mass matrix elements b-Qs and S-Q3. 
Note that whereas the Standard Model predicts a branching fraction of 6 ^ 57, including 
the next-to-leading-order correction, which is still allowed by the experimental data, future 
reduction in the experimental error with the same central value may be a potential signal 
for new vector quarks. 

In conclusion, there may still be a hint of new physics in the current precision measure- 
ments of Rfy and v4^^. If it is due to the mixing of b with heavy quarks, the only viable 
model[0 is to add a heavy vector doublet of quarks with the unconventional charges —1/3 
and —4/3. Two other models are eliminated because they require large mixings, which in 
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turn generate large shifts in ei and 63, and are thus in disagreement with present precision 
data. 
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FIGURE CAPTION 



Fig. 1. Plot of Rb versus v4^^. The solid lines are the experimental ranges. The dashed 
(dotted) line is obtained by varying gj^ {g]^) holding gj^ {g\) fixed in the Standard Model. 
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